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Abstract 

We compute the renormaUsed four-fermion operator (9^'5=2 using a non-perturbative 
method recently introduced for determining the renormalisation constants of generic 
lattice composite operators. Because of the presence of the Wilson term, O^^^^ mixes 
with operators of different chiralities. A projection method to determine the mixing 
coefhcients is implemented. The numerical results for the renormalisation constants 
have been obtained from a simulation performed using the SW-Clover quark action, 
on a 16^ X 32 lattice, at (3 = 6.0. We show that the use of the constants determined 
non-perturbatively improves the chiral behaviour of the lattice kaon matrix element 
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1 Introduction 



Renormalisation of lattice operators is a necessary step for obtaining physical results from 
numerical simulations. In this paper, we apply the general method introduced in to the 
four-fermion operator[| 

0^^=2 = (.-7^ci)(^-7,M, (1) 
which appears in the weak effective Hamiltonian relevant for K^-K^ mixing 

n^i=^ = C(M^//i)0^^=2(/i) , (2) 

where 7^ = |7^(1 — 75), 0^^^'^{ji) is the renormalised operator, C{Mw/ jj) is the corre- 
sponding Wilson coefficient and /i the renormalisation scale. The K^-K^ matrix element of 
qAs=2^^^ defines the so-called kaon 5-parameter 

{K'\0^'=\^)\K') = lflmlBM. (3) 

The uncertainty in the value of this matrix element restricts the precision with which the 
CKM matrix elements p and rj (in the Wolfenstein parametrisation) can be determined from 
experimental measurements. It is therefore of considerable importance to determine this 
matrix element using lattice simulations. 

In the continuum, chiral symmetry implies that the kaon matrix element of O^^^^ van- 
ishes in the chiral limit ^ 

{K%q)\0^'=\^,)\K%p)) = ^{p . q) + O {{p ■ q)') . (4) 

On the lattice however, in simulations based on Wilson's formulation of the fermion action 
(such as the standard Wilson action or the SW-Clover action ^), the presence of chiral 
symmetry breaking terms leads to the mixing of O^"^^^ with operators of different chirality 
[^]-||T0|, and the matrix element of O^^^^ is different from zero at p - g = |Tl[]-[0. For this 



reason, it is possible to define a renormalised operator with definite chiral properties only in 
the continuum limit, i.e. when a 0. At finite a, one can improve the chiral behaviour of 
the matrix element of O^^^"^, by subtracting a suitable set of dimension six operators. The 
mixing coefficients have so far been computed only in one-loop perturbation theory 
[ p!0[| . In this way, the systematic error in the value of the matrix element determined on the 
lattice is of 0{a'^). In addition, as a consequence of the finiteness of the lattice spacing, there 
are errors of 0{a). Following Symanzik's proposal, one can reduce these discretization errors 
from 0(a) to ©(a^a) by using the tree-level "improved" SW-Clover lattice quark action 
[^, |TB[. Using this action, the improvement has been shown to be effective for two-fermion 



operators, at values of /3 currently used in numerical simulations [0]-p2|. It remains true 
however, that ignorance of higher-order perturbative corrections to the mixing coefficients 
can distort the chiral behaviour of the operator and hence induce a large systematic error in 
the determination of Bk- The use of a non-perturbative approach to the determination of 
the renormalisation constants, should reduce this systematic effect. 



^We use the Euclidean metric throughout this paper. 
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In the following, we will define a renormalised operator 0^^^^{fi), obtained by applying 
the non-perturbative method of ref. [||] to the computation of the mixing coefficients and of 
the overall renormalisation constant. In order to reduce the discretisation errors, including 
those induced by the mixing with higher dimensional operators, it is necessary to use an 
improved fermion action and operators. In the computations described below we have used 
the improved SW-Clover action and the "improved-improved" operators introduced in ref. 
[|T0| , pO| . We monitor the effects of the non-perturbative determination of the mixing coef- 
ficients by comparing the chiral behaviour of the matrix element {K^{q)\0^^^'^\K'^{p)) jg^^|- 
computed by using the operator renormalised with standard or boosted perturbation 
theory to the matrix element of the operator renormalised non-perturbatively. In particular, 
by parametrizing the matrix element near the chiral limit in the standard way [|12]-[|15|, 



^^0|oA5=2|^o^ 1^^^ = a + f3mj, + ^{p ■ q) + 5m], + eml{p . q) + ({p . q)^ + . . . , (5) 

we investigate the differences in the values of a and /?, obtained by fitting the dependence of 
the matrix element on the kaon masses and momenta. Since in the continuum a and f3 are 
absent, cf. eq. (^, we consider a reduction of their values as a measure of the improvement in 
the chiral behaviour and in the accuracy of the determination of the matrix element. Using 
the data of the APE collaboration iT6|, we show that the chiral behaviour is indeed 



improved by using the non-perturbative results. 

The paper is organized as follows. In section 0, we briefiy summarize the strategy followed 
for computing the mixing coefficients and the overall renormalisation constant of the relevant 
four-fermion operator; in section |^, we illustrate the projection method used to determine 
the mixing coefficients. Although the method is applied specifically to the renormalisation of 
the operator O^^^^, it can readily be generalised to other sets of operators which mix under 
renormalisation. In section ^ we give some information about the perturbative evaluation 
of the renormalisation constants on the lattice; in section ^ we present the details of the 
numerical simulation and discuss our results and, finally, we present our conclusions in 
section |^. 



2 The non-perturbative method for four-fermion op- 
erators 

The renormalisation method proposed in completely avoids the use of lattice perturbation 
theory and allows for a non-perturbative determination of the renormalisation constants of 
any composite operator in a renormalisation scheme which is independent of the method 
used to regulate the ultra-violet divergences. In particular, the renormalised operators are 
independent of the fact that we start from bare operators in lattice QCD. To stress this 
point further we will refer to the renormalisation scheme defined below for 0^^^^ as the 
RI (Regularization Independent) scheme ||2^0. Non-perturbative renormalisation conditions 
are imposed directly on quark Green functions with off-shell external states in a fixed gauge, 

^ Although of course such a name could be applied equally well to many other schemes. 
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for example the Landau gauge. The method is expected to work in all cases where it is 
possible to fix the virtuality of the external states = /i^ so as to satisfy the condition 
Aqcd <^ Ai "C I /a. The condition ^ Aqcd is necessary because one has to match 
perturbatively the effective Hamiltonian, expressed in terms of operators renormalised at 
the scale /i, to the full theory. This condition is common to all approaches currently used. 
The condition /i ^ 1/a is a requirement common to all lattice methods and is due to the 
presence of 0{a) {0{asa)) effects in the operator matrix elements. The existence of the 
"window" Aqcd /i -C 1/a depends on the value of the bare lattice coupling /3 at which 
the numerical calculations are performed. We refer the reader to ref. [||] for a more detailed 
discussion on this point. 

In the following, in order to use a more transparent notation in the formulae, we will 
consider the operator 

0+ = ^[(V'i7^V^2)(^37>4) + {2^ 4)], (6) 

with four distinct quark flavours (/ = 1,2,3,4) instead of the operator O^^^^. 0+ and 
0'^^='^ have the same renormalisation properties. 

The discretization of the quark action d la Wilson, induces a mixing of the operator 
with operators of a different chirality which, in the language of refs. pi]-[El], correspond 



to the so-called "effervescent" ("evanescent") operators. The mixing, being a consequence 
of the regularization procedure, is not limited to the lattice case, but is present also in 
continuum regularizations. The effervescent operators must be subtracted from the bare one 
by a suitable renormalisation procedure |. 

In the lattice case, CPS symmetry flxes the basis of operators that may appear in per- 
turbation theory p8[| : 

Of = -Y^Pl^2)(^3^4)-(^l75V^2)(V^375V^4) + (2^4)], (7) 
0+^ = ^[(^i7m^2)(V^37mV^4) - (^17^.75^2) (^37m75V^4) 

+ (2 -'4)], (8) 

Of^ = ^^ [(V^1^2)(V^3V^4) + (V^175V^2)(V^375^4) 

+ (^lC^/..^2)(^3t^M-^4) + (2^4)], (9) 

where A^^c = 3 denotes the number of colours. 

We renormalise the operator 0+ by introducing the subtracted operator 

OM = Z+Ol = Z^{0+ + ZiOf + Z,Ol^ + ZsOf ^) , (10) 

where Z+^1^2,3 = ^+,i,2,3(Ato, 5'o('^)) ^^"^ bare lattice coupling is given by the relation 
(3 = 6/(7o(^)- The mixing constants Z,, i = 1, . . . ,3, are determined by means of projection 

^ Using dimensional regularization, the one-loop mixing with the "effervescent" operators is cancelled by 
the minimal subtraction of the pole in 1/e. 



+ ' 
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operators that will be defined in section |^. Their values are fixed by the requirement that, 
up to terms of 0{asa), renormalises multiplicatively. is logarithmically divergent as 
a — >■ 0, and this divergence is removed by imposing a renormalisation condition on which 
defines the overall renormalisation constant Z+{fia, g^la)), 

Z+{fia, gl{a))Z^^{fia, gl{a))Tl{pa)\p2^f,2 = 1, (11) 

where F^pa) is obtained by projecting a suitable amputated Green function of the operator 
0\. on the Dirac structure 7^ (see eq. ( p3| ) in section ^ and refs. is the quark 

field renormalisation constant to be defined below (eq. (|33D of section |^ and ref. In eq. 
(|TT]), = denotes the momentum of the external quark states. We have chosen equal 
momenta for all four external quark legs, because this is the simplest choice which regulates 
the infrared divergences ||2^ . 



The renormalised operator in eq. (0), calculated in the RI scheme, depends both on the 
gauge and on the external states. The Wilson coefficient must also be calculated in the same 
gauge and with the same external states in order to obtain the physical operators which 
are independent of both0. The next-to-leading order calculation of the Wilson coefficient 
relevant for the operator (|^), in the Landau gauge and with equal external momenta, can be 
found in ref. p4 . 



3 Determination of the mixing constants 

In this section, we define the four-point amputated Green functions and introduce the pro- 
jectors that have been used to determine the mixing constants. 

Since the non-perturbative renormalisation conditions are imposed on quark states, the 
Green functions of a four-fermion operator will depend on four coordinates. Denoting by 
Xi,X3 and 0:2, 2:4 the coordinates of the outgoing and incoming quarks, the Green functions 
corresponding to the insertion of the operators (^-(0) can be written as linear combinations 
of Green functions of the form 

Gr'^{xi,X2,X3,Xi) = (^/'i(xi)V52(a;2)C'r-(0)^/'3(x3)^4(x4)) , (12) 

where (■ ■ ■) denotes the vacuum expectation value, i.e. the average over the gauge-field 
configurations. The generic four-fermion operator Opa is given by 

Or.(o) = Cr« [v^i(o)r>2(o)v^3(o)r>4(o) + v^i(o)r>4(o)^3(o)r>2(o)] , (is) 

where denotes a Dirac matrix, and Cr» is a constant associated with F". The index a can 
be either single-valued (if F"^ = I or 75) or be summed over a range of values (if F*^ = 7^, 
7/i75 or a^i, a sum over repeated Lorentz indices is implied). 

The Fourier transform of the non-amputated Green function (|l^), at equal external 
momenta p, has the form 



AB ^-rat^\CD\ lTa(^\AD f-^Ta(^\CB\ 



*This is true up to higher order continuum perturbative corrections and lattice systematic errors. 
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where 

r%p)^^ = S{p\0)^Xp<r{j,S{p\0)h,)^^ . (15) 
In eqs. (0) and (15) the upper-case Roman superscripts denote colour labels and the lower 



case Greek subscripts denote spinor labels. 5(^10) is defined by 

S{p\0) = J d^xS{x\0)e-'P-'', (16) 

where S{x\0) is the quark propagator computed on a single gauge-field configuration (cf. 
section 4 of 0), and is therefore not translationally invariant. It satisfies the relation 

^(x|0) = 75^^(0|x)75. (17) 



The amputated Green function can be obtained from eq. ([I^ 



AAp)^7f = 5-'(p)?«^5-'(p)?;^Gr.(j>)i^5^5-^(p)?i^5-^(p)fJ , (18) 

where S{p) is the Fourier transform of the translationally-invariant quark propagator, i.e. 
the Fourier transform of S{x\0), averaged over the gauge-field configurations. 

As mentioned above, the renormalisation procedure necessary to determine the mixing 
constants consists in defining suitable projectors on the amputated Green functions of the 
operators (||)-(|D. To this end let us introduce a more convenient notation. Let us denote 
by Oi, i = 0, . . . , 3, respectively, the operators 0+, O^^, OX/^^ O^^^. Then, the projectors 
Pi, z = 0, . . . , 3, are defined by the condition 

TrPiAf = %, ^,j = 0,...,3, (19) 

where Af'\ i = 0, . . . ,3, are the amputated Green functions, at tree level, of the operators 
Oi, and the trace is understood over colour and spin (as defined below). The renormalisation 
scheme depends on the precise definition of the projection operators and we now define our 
procedure in detail. 

For each Dirac Matrix r**, we define the projector Pp6 by 

Tr Pr^Ar.(p) = (F^^ ® r^^,)Ar»(p)Jf/7 , (20) 

where the index b is either fixed or corresponds to a sum over repeated indices. In the free 
theory, the amputated Green function reduces to 

A?i(p)Jfp^J = CMS'^'S^'^'iV;^ ® F^,,,) - 5^^'5^'^(F^,, ® F;,J]. (21) 

and the result of the projection defined in ( |2I1| ) is: 

Tr ProA(°i(p) = Cr4N^{TT F"F'')(Tr F'^F^) - N,{Tt r^F^F^F'')]. (22) 
^It is only the traces (EG) which are required for the determination of the subtraction constants. 
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The projectors corresponding to 0+ and to the operators Oi defined in eqs.(0)-(||) are as 
follows: 



Po = ^7777^^P7?> (23) 



1 

We(iV, + 1) 
Pi = ../^.o. (Pl-P.. 



2(1 - iV2) 

1 

4(1 - iV2) 

2(l-iV2)(iV2 + iV,-l)^^I " 



+ ^(P>-P>7a), (24) 



c J \ c 

+ 4(1 - NI){n] + N,-l) ~ ^^"^-^^ ^^^^ 
Ps = g(^,^_^) (Pl + P7s+P..J, (26) 

where 7^ = §7^(1 — 75). Note that the projectors Pq and P3, eqs. ( ^ and (|26D , have the 
same Dirac structure as the operators Oq and O3, eqs. (|^) and (^. This is due to the Fierz 
rearrangement properties of these operators. 

It is possible to determine the mixing coefficients Zi by using the projectors (P5|)-(PBD 
defined in the free field case. Let us introduce the matrix D defined by 

K = Y^ D,,Af , (27) 

j=0 

where the elements Dik are determined non-perturbatively by the projections 

Afc = Tr PfcA,, t,k = 0,...,3. (28) 

The mixing constants Zi are then fixed by the condition that the subtracted operator is 
proportional to the bare free operator 



Tr PfcA; = (^Dok + E ^^Afe^ =0, A; = 1, 2, 3 , 



(29) 



Equation (|29|) yields three conditions corresponding to a linear non-homogeneous system in 
the three unknowns Zi. Defining the reduced 3x3 matrix D as 

Dik = Dki, t,k= 1,2,3, (30) 

the solutions of this linear system are given by 

Z, = -j2{D)ij}Dok, t = 1,2,3. (31) 

k=l 
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The overall renormalisation constant Z+ is then determined by the condition (pJ]), using 

riipa) = Tr PoA; = (^D^o + E ^«Ao) • (32) 



In eq. (pi]), the renormalisation constant is defined by the relation 



Z^{^a,gl{a)) = ^Tr (AyiT/.) 



X 



(33) 



where = ipjfiip is the local vector current, and ZyL its renormalisation constant which can 



be determined with high accuracy, by using the vector current Ward identities [^, ^ ^ 
AyL is defined as 

AyL{p) = S{p)-'GyL{p)S{p)-' , (34) 

where GyL{p) is the non-amputated two-point Green function of the local vector current, 

GyL{p) = (rV(p)), cf. eq. (0). There are various equivalent ways to define Z^, but (|33| ) is 
the most natural from a non-perturbative point of view. For a more thorough discussion on 
the determination of Z^, we refer the reader to sections 2 and 4 of ref. [Q]. 



4 Lattice perturbation theory 

We have also calculated the renormalisation constants Z^,Zi,Z2 and Z^ in one-loop per- 
turbation theory, in order to be able to compare the results with those obtained non- 
perturbatively. Since the non-perturbative renormalisation condition depends on the gauge 
and on the external states, the perturbative calculation must be done in the Landau gauge 
and at equal external momenta. This calculation is an extension of those of refs. ^ [To| . 

Starting from a bare lattice operator 0(a), the one-loop vertex function Tplpa) is ob- 
tained by tracing the amputated Green function (but with wave function effects included) 
with a suitable projector. The generic expression of TQ^pa), calculated between states of 
momentum p and in a fixed gauge A, is 



T'oipa) = [l + ^ (7olog(l//a2) +rJ^"(A,p)) 



(35) 



7o is the anomalous dimension, which at one-loop order is independent of the gauge, the 
external states and the regularizationQ. The finite coefficient rJ^^^(A,p), on the other hand, 
does depend on the gauge, the regularization and the external states. The momentum label 
p appearing as argument of rJ^^^(A,p) indicates that the result is a dimensionless function 
of the external states. 



^ A denotes a generic covariant gauge: A = corresponds to the Landau gauge, A = 1 correspond to the 
Feynman gauge. 
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In the continuum, in any renormalisation scheme based on dimensional regularization 
(DR=NDR, HV or DRED), the vertex function between states of momentum p and in a 
generic gauge, is 



= [i + ^ (70 log(MVp') + rg^(A,p) 



(36) 



where /z is the DR renormahsation scale. Thus, the one-loop relation between the operators 
in the continuum and on the lattice is 



0(/i) 

where 



l + ^(7olog(/.V) + A°^- 



0(a), (37) 



^DR- latt ^ ^DR^^^ _ ^ Jatt ^3g) 
is independent of both A and p. From Aq^~ \cLtt^ calculate the lattice constant 



|^^^(A,p), in any gauge and at any external momenta, from the corresponding constant 
in the continuum, rQ^{X,p). In order to compare the perturbative result with the non- 
perturbative determination, we need rJ^^^(A = 0,p) in the Landau gauge and with non-zero 
but equal external momenta. From (pH]), we immediately obtain 



^ latt 
o 



(A = 0,p)= r°^(A = 0,p) - Ag^- 1^". (39) 



Since rj^^^ must be independent of the continuum regularization used in the intermediate 
steps, a check of the correctness of the calculation is given by 

r^^^^iX = 0,p) - Ag^™- latt = rS°^(A = 0,p) - Ag^^" 1^", (40) 
which is equivalent to 

r^^^^iX = 0,p)- rg°^(A = 0,p)= rg^^°(A = l,p') - rr^(A = l,p'). (41) 

The one-loop contribution to the renormalisation constant and to the mixing coeffi- 
cients Zj's have been calculated in |P, 10|, by comparing the lattice and the DRED scheme. 
In the notation of these authors 

Z+ = 1 + F+ = Ag^™-l^** = -10.9, 

An ^ 

Zi = Z2 = Z3 = ^F*, F* = 19.4. (42) 

In the DRED scheme, for a generic gauge A and by taking the momenta of the external legs 
to be equal, we find 

rgf °(A,p) = A(-7/3 + 8/31og(2)) - 5/3 + 81og(2) , (43) 
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whilst the mixing with the "effervescent" operators is cancelled by the minimal subtraction 
procedure. Thus, the perturbative expressions of the lattice renormalisation constants in the 
RI scheme are given byj] 

Zf = l-f^(-7o,log(/iV)+rgf°(A = 0,p)-F+), 



= ^F*, (44) 

with 7o+ = —2 and rQ^™(A,p) given in eq. 

We have not evaluated the renormalisation constants using Discrete Perturbation Theory 
(DPT), i.e. by summing only over the discrete values of momenta allowed on our finite lattice 
(this was done for the two-quark operators in ref. |1[]). We have only evaluated the constants 
using standard lattice perturbation theory, in which finite lattice size effects are neglected. 

In order to estimate the values of the renormalisation constants we have used the fol- 
lowing "boosted" coupling constant |23 



ar = TTTfTTrr^s^^" ^ 1-68 (at P = 6.0) . (45) 

(3 ir Up) 

as our expansion parameter and refer to the result as corresponding to Boosted Perturbation 
Theory (BPT). We also present the values obtained using the bare coupling a}^^^, and refer 
to these results as coming from standard perturbation theory (SPT). 



5 Numerical results 

In this section, we give the numerical results of our calculation. The simulation has been 
performed by generating 36 independent gluon-field configurations, on a 16^ x 32 lattice, 
at (3 = 6.0. The errors have been obtained with the jacknife method, by decimating three 
configurations at a time. The SW-Clover quark propagators have been computed at a single 
value of the quark mass (m^a ~ 0.07), corresponding to the hopping parameter k = 0.1425. 
The quark Green functions have been computed in the lattice Landau gauge, defined by 
minimizing the functional 



Tr 



(46) 



Possible effects from Gribov copies have not been studied. For more details, see ref. |^ . 

In fig. |l], the renormalisation constants, obtained by using the prescription described 
in sections ^ and |^, are given as a function of the renormalisation scale /i^a^. We hope to 
find an interval of values of /i^a^, large enough to avoid significant non-perturbative effects 
and small enough to avoid large discretization errors. In ref. the existence of such 



^ The scale fi in this formula denotes the renormalization scale at which rQ{pa) is renormalized, i.e. 
+^ o( 
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Figure 1: Non-perturbative renormalisation constants of the operator 0+ as a function of 
/i^a^; (a) the overall renormalisation constant Z+; (h)-(d) the mixing coefficients Zj, i — 
1, 2, 3. We also report the perturbative evaluation: the dashed curve is from BPT, while the 
solid curve is from SPT. 
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28 + 05 


0.96 


0.84 ±0.03 


0.14 ±0.07 


0.30 ±0.02 


0.24 ±0.04 


1.27 


0.80 ±0.04 


0.17 ±0.05 


0.29 ±0.02 


0.21 ±0.03 


1.54 


0.82 ±0.02 


0.19 ±0.04 


0.27 ±0.02 


0.16 ±0.03 


1.89 


0.83 ±0.03 


0.22 ±0.05 


0.30 ±0.02 


0.18 ±0.03 


2.47 


0.85 ±0.02 


0.22 ±0.06 


0.33 ±0.02 


0.23 ±0.03 


SPT 


0.91 


0.12 


0.12 


0.12 


BPT 


0.84 


0.21 


0.21 


0.21 



Table 1: Values of Z^ and Z^ f« = 1,2,3 j jor several renormalisation scales ^^a?. We also 
give the results obtained at ^"^a? = 1, by using "standard" perturbation theory (SPT) and 
"boosted" perturbation theory (BPT), using an effective coupling = 1.68 a^'^". 



a "window" in /i^a^ was investigated by comparing the renormalisation constants of two- 
quark operators computed in perturbation theory with the corresponding non-perturbative 
determinations on quark states, and with the results obtained by using the Ward identity 



method |21|, |2^, |2^. In most of the cases a range of acceptable values was found in the 
interval 0.8-0.9 < /i^a^ < 1.5-2.0. At smaller values of /x^a^, in particular in the case of 
the axial current and of the pseudoscalar density (probably because of the presence of a 
pseudo-Goldstone boson contribution), the non-perturbative corrections were found to be 
large. For this reason, it is difficult to determine the renormalisation constant of the axial 
current in this way. Only at values of /x^a^ larger than 1.5-2.0, a surprisingly large value in our 
opinion, did discretization errors become clearly visible. They were signalled by the fact that 
the renormalisation constants computed at the same values of yU^a^, but with inequivalent 
components of the momentum p (e.g. p = 27r/16a(4, 4, 0, 2) and p = 27r/16a(0, 0, 0, 6)) were 
found to be different [0. This was interpreted as a signal of the breaking of the Lorentz 
symmetry due to lattice artefacts, see also ref. pll . 

For the four-fermion operators considered in this paper, we do not have the possibility 
of checking the results for the renormalisation constants by the use of Ward identities, but 
it appears that a similar situation may also occur in this case. In the region of momenta 
/i^a^ > 0.96, the renormalisation constants are determined with a relatively small error (the 
worse case being the error of Zi which is about 50% at /i^a^ = 0.96) and the dependence 
on the scale is relatively weak, as can be seen from table and fig. |l]. We notice that the 
values of Z^ and Z2 have small statistical errors even at scales smaller than 0.96 (~ 10% 
in the worst case), that Z^ has relative errors in the range of about 15-20%, and that 
Zi has the largest relative error at all the scales considered in table |l|. As for the scale 
dependence, Z2 is quite stable as a function of /x^a^, while both Zi, which suffers from the 
largest statistical uncertainty, and Z^ do not exhibit a very clear plateau, as can be seen 
in fig. |I[ In particular the value of Zi seems to increase with the scale. With the present 
statistics, we cannot determine whether the variation of Zi and Z3 with /i^a^ is real or due to 
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21 n 7^ 


69n2i 


0.96 


0.017(13) 


0.21(17) 


0.70(12) 


1.27 


0.015(13) 


0.21(16) 


0.66(11) 


1.54 


0.018(13) 


0.22(16) 


0.67(12) 


1.89 


0.023(13) 


0.22(16) 


0.69(12) 


2.47 


0.022(14) 


0.23(17) 


0.72(12) 


SPT 


-0.067(12) 


0.17(15) 


0.62(11) 


BPT 


-0.054(12) 


0.17(15) 


0.62(11) 



Table 2: Values of the coefficients a, (3 and 7 obtained from a linear fit of (0+) = 
{K^\0^^^'^\K^) iditl {P^Y ■ The results refer to operators renormalised at different scales 
/i^a^. Values of the same parameters in standard perturbation theory and in boosted pertur- 
bation theory are also reported. 



statistical fluctuations. Fortunately, as we will see below, the largest correction to the chiral 
behaviour comes from the operator O2, corresponding to Z2, which is very well determined. 
Hence the chiral behaviour of the operator 0+ is stable with respect to the uncertainties 
above. 

In order to investigate the effects of the non-perturbative corrections, we have combined 
our results with the computation of the lattice matrix elements of the four-fermion operators 
(ll)-(ID performed in ref. [l^, where a more detailed discussion of the numerical 



aspects can be found. Here we limit ourselves to a qualitative discussion of the results. 

In fig. I, we show the chiral behaviour of (0+) = {K^\0^^=^\K^) {P^Y with 
the meson at rest as a function of X = S/^f^M^/iP^)^. {P^)^ = \{0\s-f5d\K°)\^ is the 
squared matrix element of the pseudoscalar density between the meson and the vacuum. 
The variables (0+) and X are particularly convenient since they can be obtained from 
suitable two- and three-point correlation functions without any fitting procedure [1^. 

An analysis of the different contributions to the matrix element of the renormalised op- 
erator shows that the largest correction comes from the operator O2, eq. (|), whose constant 
Z2 is well determined. This contribution is much larger than that coming from Oi (and 
larger than that from O3) 0. Thus, the uncertainty in Zi (and partly Z3) has no significant 
consequences for the value of the kaon matrix element (|^). Indeed, as can be seen from table 
0, in passing from yU^a^ = 0.66 to /i^a^ = 2.47, the central value of Zi increases by a factor 
of 3, but there are no large variations in the values of a, (3 and 7, eq. (^, i.e. the "physical" 
results depend rather weakly on the scale. The use of the non-perturbative renormalisation 
constants leaves the values of (3 and 7 almost unchanged compared to those obtained by us- 
ing the constants computed in one-loop perturbation theory: (3 ~ 0.2, and within the errors 
is compatible with zero in both cases, and 7 is about 15% larger in the non-perturbative 

^Notice that Oi has the smallest colour factor. 



12 



0.20 



0.15 



0.10 



0.05 



0.00 



-0.05 



T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



<0^> 



-E3- 



<> 



[] 
<> 



o SPT □ BPT o NP 



_Q I I I I I I I I I I I I I I I I I I I I 

0.05 0.1 0.15 0.2 



X 

Figure 2: Chiral behaviour o/(0_|_) as a function of X (see text). We give the matrix elements 
of the operator renormalised in standard perturbation theory (O), boosted perturbation theory 
(D) and non-perturbatively (o). 



case. In contrast, a changes sign and its absolute value is reduced by about a factor of 
three in the non-perturbative case, and becomes compatible with zero. This happens for 
any choice of /z^a^ between 0.46 and 2.47, cf. tab. |l]. The variation of the Zj's in the interval 
of /x^a^ considered in tables |I] and |^ is representative of the variation allowed by the statis- 
tical errors. Since the Zj's and the matrix elements of the four-fermion operators have been 
computed on different sets of configurations, this is the best test of the stability of the results 
which can be done at present. A more definite conclusion will be reached by computing the 
renormalisation constants and the matrix elements on the same set of configurations. 

Since a should vanish in the continuum limit, we conclude that the use of the non- 
perturbative renormalisation constants improves the chiral behaviour for a large range of 
values of the renormalisation scale. This is also illustrated in fig. ^ obtained for 



/i a = 0.96. More details will be given in ref. |17 . 
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6 Conclusions 



We have applied the non-perturbative renormahsation method proposed in ref. to the 
AS* = 2 operator given in eq. (|1]). Since, on the lattice, this operator mixes with other 
dimension-six operators of different chirality, we have illustrated a projection method for 
the determination of the mixing coefficients. The overall renormahsation constant of the 
subtracted operator has then been obtained as in the case of any other multiplicatively 
renormalisable operator. 

In this exploratory study we have computed the subtraction constants with limited 
statistical precision (36 configurations on a 16^ x 32 lattice at (3 = 6.0 using the improved 
SW-Clover action). The results in fig.|I] and table |I| are very encouraging, and motivate us to 
repeat the calculation with larger statistics and at different values of (3 {(3 = 6.2 and 6.4), and 
to extend it to the operators relevant for A/ =1/2 transitions and to the penguin operators 
which control CP-violation in kaon systems. Even with our limited statistical precision, 
our results indicate that the chiral behaviour of the K^-K^ matrix element of 0^^^^{fj,) is 
improved significantly by the use of the subtraction constants which were determined non- 
perturbatively. This supports our view that, by combining the improvement of the action d 
la Symanzik, which reduces 0{a) effects, with the non-perturbative method of ref. [|1|, which 
reduces higher-order effects in the mixing coefficients, it is possible to achieve an accurate 
determination of the physical weak amplitudes using Wilson-like fermions. 
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